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Abstract. We evaluate the validity of leading models of the Galactic magnetic field for 
predicting UHECR deflections from Cen A. The Jansson-Farrar 2012 GMF model (JF12), 
which includes striated and random components as well as an out-of-plane contribution to 
the regular field not considered in other models, gives by far the best fit globally to all-sky 
y—i data including the WMAP7 22 GHz synchrotron emission maps for Q,U and ~ 40,000 

extragalactic Rotation Measures (RMs). Here we test the models specifically in the Cen A 
region, using 160 well-measured RMs and the Polarized Intensity from WMAP, nearby but 
outside the Cen A radio lobes. The JF12 model predictions are in excellent agreement with 
the observations, justifying confidence in its predictions for deflections of UHECRs from 
Cen A. We find that up to six of the 69 Auger events above 55 EeV are consistent with 
originating in Cen A and being deflected < 18°; in this case three are protons and three 
have Z = 2 — 4. Others of the 13 events within 18° must have another origin. In order for 
a random extragalactic magnetic field between Cen A and the Milky Way to appreciably 
alter these conclusions, its strength would have to be > 80 nG - far larger than normally 
imagined. 
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1 Introduction 

The astrophysical origin of the highest energy cosmic rays remains one of the major open 
questions in astrophysics. Cosmic rays with energies above about 55 x 10 18 eV (55 EeV) - 
commonly called UHECRs - lose energy through interactions with the Cosmic Microwave 
Background radiation, implying that their sources must be relatively nearby (closer than 
~ 200 Mpc[l, 2]). However, very few classes of astrophysical systems seem capable of accel- 
erating particles to such high energies and there are few candidate sources within the required 
distance. The closest plausible UHECR source is the nearby radio galaxy Centaurus A[3-5]. 
Tantalizingly, the highest energy Pierre Auger Observatory events show a distinct excess 
within 18° of Cen A: 13 events observed, while 3.2 events would be expected for 69 events 
coming from an isotropic distribution [6]. This excess has re-ignited interest in the possibility 
that Cen A is the source of most of these events[7-13]. However, Cen A is located in the 
"Supergalactic Plane" - a nearby overdensity of galaxies - so the UHECR excess could be 
produced by multiple sources in the general region, or possibly simply be due to some fo- 
cussing effect of the GMF or EGMF. Fig. 1 shows the arrival directions of the 69 published 
events, with Cen A and the supergalactic plane indicated as well. 

Since UHECRs are charged protons or nuclei[14], their arrival directions do not point ex- 
actly toward the source of their emission. Magnetic deflection is proportional to the UHECR's 
electric charge, which is undetermined. Some evidence, e.g., correlations with extragalactic 
objects such as AGNs [6, 15], favors UHECRs being mostly protons, while other evidence 
suggests a mixed or heavy composition [16]. Therefore, use of a reliable model of the Galac- 
tic magnetic field including both the large scale and random components, is necessary to 
obtain trustworthy predictions for UHECR source locations. This is not a pedantic matter, 
as demonstrated in Fig. 2 which shows the striking variation in predicted arrival directions 
using six currently-used models for the coherent GMF, for a 60 EeV proton produced in the 
core of Cen A. 

Of the six GMF models used in Fig. 2, the new model of Jansson and Farrar [17, 18] 
(JF12 below) gives by far the best global fit to the RM and polarized synchrotron data [17], 
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Figure 1. The 69 published Auger events above 55 EeV[6] in Galactic coordinates, with Galactic 
longitude zero at the center and increasing to the left. The giant radio lobes of Cen A are marked in 
red and the super-Galactic plane is shown as a grey line. The blue contour is an 18° circle around 
Cen A. The Auger exposure is shown in light blue; it is largest at the South Celestial Pole, about 50° 
south of Cen A. 

with a Xdof P er degree of freedom of 1.096 for the 6605 observables (pixels of RM, Q and 
U). The next best model is that of Sun et al.[19], whose parameters were updated in [20] 
(SR10). This is the most comprehensive attempt prior to JF12 to model the coherent GMF 
using constraints from both RM and synchrotron emission data. It does not allow for as 
general a functional form as JF12, and in particular does not include the out-of-plane field or 
possible striated random fields. With the parameters given in ref. [20], its Xdof f° r ^ ne 6605 
observables is 1.672. To test if the functional form adopted in SR10 is as good as JF12, we 
used the SR10 form and re-optimized its parameters to fit the ensemble of the data using the 
JF12 MCMC; the resulting fit has vj f = 1-325, indicating that the out-of-plane and striated 
field components of JF12 are signficant improvements in the model. Although more recent, 
the models proposed by Pshirkov et al.[21] (P+ll) have a less general form and are less well 
constrained, as these authors used only RMs and not synchrotron emission. They are based 
on the SR10 and Prouza-Smida[22] (PS03) models. Being unable to disambiguate the large 
scale geometry, P+ll offers benchmark BSS and ASS versions. When fitting the complete 
set of 6605 observables, these give Xdof = 2-663, 4.971, respectively; with our re-optimization 
of their parameter values these become Xdof = 1-452, 1.591. We did not measure the x^of 
of PS03 because the P+ll models are a generalization of it, and those give poor fits. The 
ASS and BSS models of Stanev[23] are classics, developed to illustrate the impact of different 
field geometries more than to provide a detailed model for the field; they fare even worse in 
a global fit. Clearly, studies of the deflection of UHECRs such as refs [24, 25] using these old 
models, cannot be trusted to reliably predict CR deflections in the direction of Cen A. 

In Sec. 3 we take advantage of the exceptional RM coverage in the region surrounding 
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Figure 2. The predicted locus of arrival directions for a 60 EeV proton emitted from the nucleus 
of Cen A (white circle), for the JF12 GMF and five other popular large-scale Galactic magnetic field 
models: the ASS/BSS models by Stanev [23], the best-fit model of Sun et al. [19, 20] with a 2 ^G 
halo field, and the ASS/BSS models of Pshirkov et al. [21]. The 2 — a uncertainty region of the 
predicted arrival direction due to the uncertainty in the JF12 parameter values is indicated by the 
shaded region; no such uncertainty analysis exists for the other models. JF12 provides a model of the 
random field, but for purposes of comparing to the other models which do not provide a model of the 
random field, only deflections due to the coherent field are shown. 



Cen A from [26], to test the various GMF models in the region relevant for predicting 
deflections of UHECRs from Cen A. We find that JF12 accurately predicts the mean Faraday 
rotation measure and polarized and total synchrotron intensity in the particular direction of 
Cen A, while other models perform less-well to very-poorly. 

Finally, having confirmed the validity of the JF12 model for Cen A deflections, we use 
JF12 in Sec. 4 to determine the deflections of UHECRs through the GMF as a function of 
their energy and charge. We find that three events within 18° of Cen A could be protons 
coming from Cen A and three others can be attributed to Cen A for more general charge 
assignments. Thus we find that the distribution of the arrival directions of the excess of 
events is not compatible with their dominant source being either the Active Galaxy or the 
extended radio lobes of Cen A, unless high Z nuclei can "wrap back" to the Cen A region - 
winding up arriving from that direction after deflections greater than 2ir. Of course, in that 
case, an association with Cen A would be essentially accidental. 
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2 Centaurus A 



Centaurus A (NGC 5128) is the nearest active galaxy, and a Fanaroff-Riley Class I (FR-I) 
radio galaxy (see [27] for a review), at a distance of 3.8 Mpc [28]. The massive elliptical host 
galaxy has Galactic coordinates (I, b) = (309.5°, 19.4°). Thanks to its proximity and size, 
its enormous radio lobes combine into the largest extragalactic object on the sky, with an 
angular size of 9°x5°, corresponding to a physical size of 500x250 kpc. About 5 kpc from 
the central galaxy, jets from the accretion disk surrounding the central supermassive black 
hole expand into plumes as they plow into the ambient intergalactic medium. These plumes 
are called the inner radio lobes. Some material goes farther, creating the northern middle 
lobe, which extends to 30 kpc and lacks a southern counterpart. The giant outer radio lobes 
extend 250 kpc in projection both in the north and the south; their outline is shown in Fig. 
3. The 3D orientation of the lobes is not well-known. 

Cen A was first considered as a possible source of UHECRs by Cavallo [3]. The possi- 
bility that Cen A could in fact be the source of most cosmic rays, if turbulent extragalactic 
magnetic fields on the Cen A side of the Milky Way are near the maximum allowed value, 
was first proposed by Farrar and Piran [5]. Some of the more recent works include: [7], which 
proposes that Cen A is one of three sources that, combined, are responsible for all observed 
UHECRs; [8] which analyzes the significance of correlation between UHECRs and Cen A; [9] 
which investigates the plausibility of the giant radio lobes of Cen A being acceleration sites 
of UHECRs; [10] which considers the possibility that the radio jet at the core of Cen A is 
an accelerator; [11, 29] which considers various mechanisms that could accelerate UHECRs 
in radio galaxies such as Cen A; and [12, 13] which argue that Cen A is the source of the 
~ 10 events in the region of the sky surrounding Cen A. Ref. [30] notes that Cen A could be 
associated with at least 4 out of the - at the time - 27 published Pierre Auger UHECRs above 
57 EeV due to its large radio extent, but does not consider the deflection by any particular 
Galactic magnetic field model. 

3 Magnetic field model and predictions 

As seen in Fig. 2, different GMF models yield highly disparate predictions for UHECR 
deflections. Hence, the validity of any conclusions about UHECR deflections hinges on the 
reliability of its GMF model. The JF12 model of refs. [17, 18] should be substantially more 
realistic than earlier models. It has a more general form and it is constrained by both RMs and 
synchrotron emission, which together probe both the line-of-sight and transverse components 
of the field. The model includes a thin disk component, an extended halo component, and 
an out-of-plane component as suggested by observations of external galaxies; random and 
striated random fields are also included in the analysis. We refer the reader to [17] for details 
of the JF12 large scale GMF model, and to [18] for a description of the associated random 
field model. The model is constrained by the WMAP7 Galactic synchrotron emission map 
[31] and more than forty thousand extragalactic rotation measures, and as noted in Sec. 1 it 
reproduces the global RM and polarized and unpolarized synchrotron data well. 

However, the Galactic magnetic field is very complicated and even the JF12 global 
model with 34 parameters describing the coherent and random fields cannot be expected to 
provide a highly accurate model of the magnetic field along every line-of-sight. Therefore, 
before proceeding to UHECR deflections, we first determine the constraining observables 
along lines-of-sight relevant for UHECRs propagating from Cen A and compare them to the 
predictions of the JF12 model. 
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Figure 3. Polarized synchrotron radiation at 22 GHz (color) from WMAP7 data [31]. The published 
Auger UHECR events above 55 EeV [6] around Cen A are indicated with small gray circles; their 
energies are given in EeV. Contours from radio data [33] outline Cen A (center) and parts of the 
Galactic plane (both contours drawn at 70 K at 408 MHz). 160 extragalactic sources with line-of- 
sights outside Cen A [26] are shown along the boundary of Cen A. Filled white circles denote negative 
rotation measures (corresponding to a line-of-sight electron-density-weighted average magnetic field 
away from the observer), black squares (very few) denote positive rotation measures. The size of 
the markers is proportional to the magnitude of the rotation measure. The large white circle shows 
the region 5° to the right, used to estimate the PI in the direction of Cen A without foreground 
contamination. 

Fig. 3 shows a portion of the sky centered on Cen A, with color scale indicating the 
synchrotron polarization intensity (PI) from WMAP7. It can be seen from Fig. 3 that Cen A 
is located at the edge of a highly polarized region, part of the nearby North Polar Spur (NPS) 
or radio loop 1[32]. If the measured PI near Cen A is dominated by these local structures, 
the emission is likely coming from very nearby (< 200 pc in the case of the NPS), in which 
case this region is not a good indicator of the large scale magnetic field relevant for UHECR 
deflection. The measured PI in the direction of Cen A is likely dominated by emission from 
Cen A itself, and also cannot be used. But the area immediately to the right of Cen A 
appears mostly uncontaminated by local emission. We select the 2° radius disk shown in Fig. 
3, centered on I = 304.5°, b = +19.4°, a point 5 degrees from Cen A; this point was chosen 
to be close to the direction of Cen A yet avoid any obvious contamination due to foreground 
sources. Here we compute the average PI and I and their variance. The PI is so small that 
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the polarization angle cannot be reliably determined and we do not use it. In this way we 
estimate the polarized intensity and intensity at 22 GHz in the general direction of Cen A to 
be PI = 0.008 ± 0.006 mK and I = 0.14 ± 0.07 mK, where the uncertainties quoted are the 
standard deviation of the l°xl° subpixels. 

The rotation measures of 160 extragalactic sources (EGS) with lines-of-sight near, but 
outside, Cen A were measured using the Australian Telescope Compact Array by Feain et 
al. [26] (see Fig. 3). Nearby foreground structures do not disproportionately impact RMs 
as they do for synchrotron emission, so these RMs provide an excellent measurement of the 
RM in this direction. The average RM of these 160 new measurements is — 54radm -2 , with 
a standard deviation of 32radm~ 2 . 

The corresponding predictions of the combined JF12 coherent and random GMF models 
are PI = 0.00895 mK, I = 0.17 mK and RM = -51.092 radm" 2 . Although the JF12 model 
parameters were determined using all available RMs including those of Feian et al.[26], the 
fitting procedure reduces the 160 RMs surrounding Cen A to just a handful of data points (4°- 
by-4° pixels) out of several thousands used to constrain the GMF parameters, so there is no 
a priori guarantee that the model prediction will agree well with the local observables. Thus 
the virtually perfect agreement between the JF12 model predictions and the observations 
along the Cen A line of sight, is a testimonial to the quality of the functional form of the 
model and the power of the global fitting approach. 

The uncertainty on the mean is 1 /y/~N times the standard deviation, where N is the num- 
ber of statistically independent measurements. The extent to which the different subpixels 
for PI and the individual RMs can be considered independent measurements depends on the 
scale length of the random fluctuations, which has not yet been established. Even with iV's of 
16 and 160, leading to observational mean values of PI = 0.008±0.0015, I = 0.14±0.025 mK 
and RM = —54 ±3, the JF12 predictions are in excellent agreement within the errors. Given 
the much worse global fit of other GMF models to the data, it is not surprising that they also 
give a poorer fit in the Cen A direction with {PI (mK), RM = (radm -2 )} predictions: SR10 
{0.00656, -25.486}, P+11BSS {0.0092, -72.411}, P+llASS {0.0096, -100.469} , Stanev 
BSS {0.00426, -56.594}, Stanev ASS {0.00426, 61.704}. 

One might think that polarized synchrotron data could be used directly to determine the 
deflections of UHECRs along a line-of-sight of interest, since both depend on the components 
of B transverse to the line of sight, but this is not the case for several reasons. Most 
importantly, the Stokes parameters Q and U are sensitive only to the orientation and strength 
of the transverse field, but not its sign, whereas a UHECR's deflection depends on the sign of 
B±: a trajectory which traverses a region with a sign reversal has a smaller net deflection than 
in the absence of the reversal. In addition, Q and U are weighted by the relativistic electron 
density, n cre , making them much less sensitive to the field in the halo than are UHECR 
deflections. (This means, of course, that when UHECR sources and charge assignments are 
known, UHECR deflections will become a very important tool for constraining the GMF.) 
Therefore it is a necessary but not sufficient condition that a GMF model gives a good fit 
to the PI data along its trajectory: a good global fit to the observables is also necessary to 
properly constrain the sign reversals. Fig. 4 shows the coherent field along the line of sight 
to Cen A in the JF12 GMF model; 9 = is angle of B± with respect to the Galactic plane, 
in a right-handed coordinate system with z pointing in the direction of motion. 
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Figure 4. The coherent field along the trajectory from Cen A, in the JF12 model. 



4 Deflections of UHECRs from Cen A 



If there were no random component to the deflections of UHECRs, the arrival directions 
of UHECRs from any individual source would fall on a single arc of width the angular 
resolution of the Auger events (~ 1°), with the separation from the source region inversely 
proportional to the rigidity of the cosmic ray {E/Z = energy /charge), in the limit of small- 
angle deflections. This distribution of arrival directions is manifestly not observed in the Cen 
A region. However, if coherent fields are negligible compared to random ones, the arrival 
directions would be smeared over all directions rather than aligned. Thus, the viability of 
the proposal that Cen A is the source of Auger's reported excess of UHECRs within 18° with 
respect to the isotropic expectation, depends on the relative importance of coherent and 
random deflections and the magnitude of the latter. Therefore it is crucial to our analysis 
that the random fields along the trajectories of interest are well constrained as has been done 
by JF12 [18]. The effect of Galactic random fields on UHECR deflections was first considered 
in [34]; with the better-constrained JF12 random field, especially in the Cen A direction, the 
uncertainty in the random deflections is reduced. 

JF12 also considered the possibility of striated fields - a field with a definite orientation 
but having a random sign and magnitude, as could be produced by stretching a random field. 
While random fields cause dispersion in all directions on the sky, striated fields aligned with 
the local regular field merely rescale the displacement due to the coherent field, similarly to 
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Figure 5. Colored regions indicate the locus within which UHECRs should be observed if the source 
is the core of Cen A; green, red, blue and ochre regions are for UHECRs of rigidity 160, 80, 40 and 
20 EeV/Z, respectively; the super-Galactic plane is shown as a grey line; nearby UHECRs and their 
energies are indicated, those with open circles are more than 18° from Cen A. The dispersion due 
to the random (turbulent) component of the magnetic field and observational uncertainties are both 
included in producing these probability distributions. Note that the actual arrival directions will fall 
within a much narrower region within this domain because the trajectories of UHECRs of similar 
energies probe approximately the same random field. The uncertainty in the mean deflection from 
the regular field, due to 1-sigma uncertainties in the GMF parameters, is shown in Fig. 2. 

UHECR energy resolution. Using RM, Q and U alone, the existence of a striated random 
field aligned with the local coherent field cannot be distinguished from a rescaling of the 
relativistic electron density, n cre , and [17] only constrained the product. The degeneracy can 
be broken by fitting the total intensity as is done in ref. [18]; we incorporate the effect of the 
striated random field along with those of the normal, isotropic random field, in the deflection 
calculations below. 

4.1 Deflections in the Galactic magnetic field 

The deflection angle is inversely proportional to the cosmic ray rigidity, for small deflections. 
Using the JF12 best-fit global GMF model, we find that the magnitude of the arrival direction 
deflection, of a UHECR from the direction of Cen A in the regular magnetic field, is 

A6 rcg = (2.3° ± 0.24°) (Z/E wo ), (4.1) 
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valid for energies for which this is a small deflection. Here, Z is the charge in units of 
the proton charge, -Eioo is the energy of the UHECR in units of 100 EeV, and 0.24° is the 
uncertainty due to GMF parameter uncertainties (the standard deviation of the separation 
in arrival directions relative to that of the best-fit-field, using GMFs for randomly chosen 
parameter sets drawn from the Markov Chain Monte Carlo probability distribution of JF12 
parameters). The CR is deflected so that it arrives from a direction closer to the Galactic 
plane and somewhat farther from the Galactic center, as shown in Fig. 5. 

To estimate the locus of uncertainty in arrival direction due to random magnetic de- 
flections we repeatedly propagate UHECRs of different rigidity through the Galaxy, dividing 
the path into domains of size A. The field in the ith domain is the sum of the JF12 coherent 
GMF, Bi, plus a random part, B raTlt x^h, and a striated-random part -Bstri,^, where h is a ran- 
domly oriented unit vector chosen to have a different random direction in each step and f\ is a 
unit vector oriented along the direction of the local coherent field with sign randomly chosen 
in each step. -Brand,* is the rms random field evaluated in the center of the ith domain using 
the JF12 random field model[18], and S s tri,i is yffi\ B m g, % \ , with (3 = 1.38[17, 18]. The domain 
size A corresponds to the maximum coherence length of the turbulent field; this is uncertain 
but is plausibly of order ~ 100 pc, a typical maximum size of supernova remnants[35, 36]. As 
one would expect, the spread in arrival directions oc vA- For UHECRs arriving from Cen A 
the mean angular separation in 1000 different realizations of the random field, between the 
arrival direction and the centroid of the arrival-direction distribution, is 

""rand — 1.3°(Z/ J B 100 )v / A^, (4.2) 

where Aioo = A/100 pc. A more sophisticated treatment of the random field realization with 
a Kolmogorov spectrum would reduce the dispersion compared to Eq. 4.2, since in that case 
power is shared over a range of scales rather than being concentrated in the largest coherence 
length which is most effective at deflecting UHECRs. 

It is important to emphasize that the dispersion quoted in Eq. 4.2 and the colored 
regions shown in Fig. 5 do not represent the typical spread of arrival directions predicted 
for UHECRs coming from Cen A, because UHECRs of a given energy from a given source 
direction and entry point into the GMF follow nearly the same trajectories and thus propagate 
in nearly the same fields for much of their trajectories, whereas the calculation above takes 
different fields for each UHECR to obtain the inclusive region of possible arrival directions, 
given our knowledge of the GMF. We can gain insight into whether different events of the 
same energy probe different coherence regions as follows. The deflection of an individual 
UHECR in a given random magnetic field can be described as a random walk with net 
deflection y/~N times the deflection 5 due to the random field in a typical coherence length: 
°"rand ~ V^V 5, where y/N is the number of independent random domains. From Fig. 4 and 
taking a 100 pc coherence length, yN ~ 12 crossing the Galaxy en route from Cen A. With 
""rand ~ 1-3° we infer 5 ~ 0.1°, so the lateral displacement relative to the trajectory in the 
coherent field alone, in crossing a typical coherence region, is ~ 0.2 pc. This is far smaller 
than the assumed 100 pc size of the coherence region, so random UHECR deflections do 
not cause trajectories for events of the same energy to diverge enough to sample different 
random magnetic field domains. More study is needed to determine whether the divergence 
of the beam from Cen A or the lateral size of the portion of the beam focussed on Earth, 
may be sufficient for different events to probe different magnetic domains. Thus Eq. (4.2) is 
an upper bound on the magnetic dispersion from random Galactic fields in the JF12 model. 
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4.2 Deflections in extragalactic magnetic fields 



The direction of arrival to the Galaxy is also smeared compared to the direction of the source, 
by deflections in the Cen A system itself and by the extragalactic magnetic field (EGMF) 
between Cen A and the Milky Way. But even very large deflections within the Cen A region 
can only move the apparent source direction within the radio lobes, translating the image 
region accordingly. Thus we are concerned with the EGMF. 

The EGMF is thought to be turbulent, with a field strength typically assumed to be 
of order 0.1 - 1 nG [37] and coherence length Aeg ~ 100 kpc. It was pointed out in [38] 
that the EGMF of the low-redshift universe is poorly constrained by the dispersions in RM 
of high redshift sources as was done in the classic work of [37], and also that [37] uses 
an unrealistically large value for the electron density and thus underestimates the true field. 
Refs. [39-42] provide an overview of recent efforts to constrain the EGMF using TeV blazers. 

If UHECRs experience many small deflections over a distance D from source to Galaxy, 
in a turbulent EGMF whose rms value and coherence length are -Beg an d Aeg; then the rms 
angle between the source and the UHECR arrival direction to the Galaxy is [43] 



assuming the lateral extent of the EGMF is > D 59. The extragalactic random deflection 
predicted in Eq. (4.3) must be added in quadrature with the Galactic random deflection 
given in Eq. (4.2) to obtain the total random deflection. In order for <5#eg ^ 12° (or just to 
be greater than the coherent Galactic deflection in Eq. 4.1), we would need 



Deflections in a variety of realizations of extragalactic fields between Cen A and the edge 
of the Galaxy were recently considered in [44] . They attempted to find choices of rms strength 
and maximum coherence length of the EGMF in the intervening region, which could produce 
the observed distribution of UHECR arrival directions within 18° of Cen A, assuming that 
Cen A is the source and ignoring GMF deflections. They were interested in the possibility 
that smaller rms field strength and larger coherence length than we employed above, could 
generate the observed pattern. However they find that except for specially selected examples, 
EGMFs do not reproduce the combination of small mean deflection and large rms deflection 
seen in the Auger events within 18° of Cen A, taking Cen A to be the source. 

4.3 Futuristic Note: RMs of Cen A sources and UHECR deflections 

In the small-angle deflection approximation one can directly relate the smearing in UHECR 
arrival directions to the dispersion in RMs accumulated in any given region in which the 
electron density is roughly constant. If the magnetic field is turbulent with a typical coherence 
length A and rms strength B, over a region of size R = N A, the rms dispersion in RM is (in 
rad m -2 ) 




(4.3) 




(4.4) 




(4.5) 
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where n e is the electron density in cm 3 . A UHECR of energy i?ioo EeV and charge Z 
propagating from the center of such a region would typically have an angular deflection 




/armor 



A 



0.52 x \ZiVAioo 




(4.6) 



where the factor y2/3 is due to the projection of the Larmor radius to the plane transverse 
to the line-of-sight. Combining Eq. (4.5) and (4.6) we obtain 



Unfortunately, in spite of the greatly improved data on RMs near Cen A [26], Eq. (4.7) 
does not provide a useful constraint on the extragalactic random field between Cen A and 
the Milky Way, because the Galaxy's much larger electron density and field strengths means 
that Faraday rotation in Galactic fields dominates the RM accumulated between Cen A and 
the Galaxy. When the GMF and electron densities are much better known, Eq. (4.7) could 
in principle become useful to eliminate the sensitivity of the random deflection prediction to 
the uncertain coherence length and field strength. 

5 Cen A as the source of proton UHECRs 

In the JF12 GMF model, a 60 EeV proton is deflected ~ 3.8° toward the Galactic plane, at an 
angle of about 45° away from the Galactic center. This prediction has an rms uncertainty of 
~ \ about the mean deflection for the best-fit GMF parameters, coming from the uncertainty 
in GMF parameters; the locus is shown in Fig. 2. Smearing due to random magnetic 
deflections is sub-dominant by at least a factor of two compared to the deflection due to the 
regular field, as can be seen in Fig. 5 which shows the maximal arrival locus of UHECRs 
of rigidities 160, 80, 40 and 20 originating in the Cen A galaxy. The Auger collaboration 
estimates a ~ 14% statistical and ~ 22% systematic uncertainty on the UHECR energies 
[45]. Adding them in quadrature produces a ~ 26% energy uncertainty for each event. Since 
deflections ~ this produces an additional uncertainty factor in the overall deflection 

ranging from 0.8 — 1.35. 

Taking a generous view of the random deflections and uncertainties in the GMF, it 
appears that up to three of the 69 published Auger UHECRs above 55 EeV could be protons 
originating from Cen A or its radio lobes. Two to three more events could be consistent with 
a Cen A origin, taking into account that the GMF model is probably not perfect, if their 
charges are Z ~ 2 — 4; events with higher charges would be deflected more than 18°. 

The Galactic magnetic field deflects UHECRs coming from Cen A or its radio lobes into 
a swath, with events of lower rigidity further from the source. Fig. 5 shows the expected 
swath of UHECRs for rigidities, E/Z = 20, 40, 80, and 160 EV if Cen A itself is the source. 
The opening angle of the region containing the UHECRs is given by o" ran d/| A# reg | < 35°. As 
long as E/Z is small enough that <7 ran d and A# reg do not become large, the numerator and 
denominator in this expression are both inversely proportional to the rigidity of the cosmic 
rays. Thus the opening angle of the swath is independent of the energy calibration, charge, 
or composition, of the UHECRs. 

The maximum opening angle of the arrival-direction-swath for UHECRs originating in 
Cen A is not only independent of uncertainties about the UHECR, it is a particularly robust 
and reliable prediction of the JF12 model. This is because the largest source of uncertainty 
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in the JF12 work is its dependence on n cre (r) which is based on GALPROP modeling[17, 
18]. But n cve {r) is a factor common to all of the observables /, Q, and U constraining 
respectively the random field and transverse coherent field, so the relative strength of the 
coherent and random deflections and hence opening angle of the swath, is quite robust. 
Details of the coherent field will change as re cre and the functional form of the field model 
are refined, resulting in possible rotation of the orientation of the arrival-direction-swath. 
However the general relationship between magnitudes of random and coherent deflections 
is likely to persist as noted, so the swath would have a roughly similar appearance, but 
possibly be pivoted about Cen A by some amount and/or the footprint be slightly stretched 
or squished. 

Hence, it is unlikely that more than 5 to 6 of the 13 Auger UHECRs within 18° in 
the published 69 event data-set can be attributed to Cen A or its radio lobes, in spite 
of uncertainties in the composition or energy calibration. We emphasize that the arrival 
direction locus shown in Fig. 5 indicates the probability distribution for arrival directions for 
the given random field: it does not represent the actual image of Cen A in UHECRs, because 
the actual GMF is a particular realization of the random ensemble used to produce Fig. 
5. The actual arrival directions should fall in a considerably narrower band, whose width 
depends on the spatial properties of the random field, especially its maximum coherence 
length, which are not yet determined. It is tantalizing that 5-6 events do fall on such a 
narrow arc within the predicted Cen A locus. 

6 Summary 

We have checked that the new JF12 model of the Galactic magnetic field [17], which gives a 
very good global fit to a large amount of Faraday Rotation Measure and synchrotron emission 
data, also gives an accurate accounting of the extensive RM data as well as the synchrotron 
data in the particular direction of Cen A. This justifies confidence in the predictions of the 
JF12 GMF model for UHECR deflections. 

Using the JF12 model of the large scale GMF[17] and the random and striated fields [18], 
we determine the locus in which protons and low-Z nuclei from Cen A should be found (Fig. 
5). Three UHECRs in the published Auger events above 55 EeV (those closest to Cen A), 
have arrival directions consistent with their being protons which originated in Cen A or its 
radio lobes. Three more events with energies 58, 78 and 64 EeV fall in the arrival locus for 
rigidity E/Z < 20 EV originating from Cen A; they could have Z = 2 — 4 and originate in Cen 
A. Events with higher charges and E < 60 EeV are deflected more than 18°. Extragalactic 
fields would have to be ~ 80 nG - far stronger than conventionally assumed - in order for 
most UHECRs within 18° of Cen A to have been produced by Cen A. 
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